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ABSTRACT 


The rich biodiversity of frogs and toads in Bangladesh is unexpected 
considering the rather simple topographic features of the country. Indeed, a 
new species, Microhyla nilphamariensis, has recently been described from 
Nilphamari District (northern region), Bangladesh. However, additional 
genetic, acoustic, and morphometric data are necessary to precisely delineate 
the newly described species as well as to measure conservation management. 
In this study we performed genetic, acoustic, and morphometric analyses. Our 
analyses showed that M. nilphamariensis was genetically divergent from its 
near congener M. ornata at 3.2% and 10.6% for the 16S rRNA and Cytb 
genes, respectively. Mean call duration of M. nilphamariensis was 0.42 + 0.01 
s (n = 8) and the call was composed of about 15.13 + 0.35 rapidly repeating 
pulses with a pulse rate of 37.9 + 0.4/s. Dominant frequency bands of M. 
nilphamariensis were much higher than those of M. ornata and M. fissipes. 
Principal component analyses showed that M. nilphamariensis differed from 
its near congeners in having shorter first and second fingers, shorter first toe, 
and longer inner and outer metatarsal tubercles relative to the snout-vent 
length. 


KEYWORDS: Mitochondrial DNA, Advertisement call, Morphometry, 
Microhyla, Bangladesh 


INTRODUCTION 


The Microhylidae is a large family of frogs and toads in order Anura 
that comprises 7.7% (572 species) of all frogs (AmphibiaWeb, 2015). 
Members of the family are distributed across Asia from the Ryukyu 
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Archipelago in Japan and China to the north, through India and Sri Lanka in 
the southwest, and through Southeast Asia to Sumatra, Borneo, Java, and Bali 
(Frost, 2015; AmphibiaWeb, 2015). Despite the wide distribution, only six 
nominal Microhyla species (Microhyla ornata, M. berdmorei, M. rubra, M. 
mukhlesuri, M. mymensinghensis, and M. nilphamariensis) are known to occur 
in Bangladesh (Kabir et al., 2009; Hasan et al., 2014a; Howlader et al., 2015). 
The accurate description of a newly described species is an Important element 
in both ecology and evolutionary biology research, in particular, in the 
development of biodiversity and conservation management strategies and 
plans (Bickford et al., 2007). 


During herpetological fieldwork in Bangladesh conducted over a 
decade, we first reported M. ornata-like specimens (IABHU 22135-22137) 
from Parbatipur in Dinajpur District (Hasan et al., 2012a). Later, as a 
continuation of our survey we only found this species from Berakhuti, Barua 
(Nilphamari District) in the northern region of Bangladesh. Therefore, it seems 
that this species may be restricted to Dinajpur and Nilphamari districts. We 
reported on similar endemicity of other microhylid frogs previously (Hasan et 
al., 2014a). 


Howlader et al. (2015) recently described M. nilphamariensis (the 
species that was identified previously as M. cf. ornata by Hasan et al., 2012a) 
from Nilphamari, based on 16S rRNA gene (165) and morphological data. 
However, additional cytochrome b gene (Cytb), acoustic, and morphometric 
data are necessary to increase our understanding of M. nilphamariensis 
because a species is the fundamental unit of biology (Mayr, 1982) and it is not 
sufficient to define biological species on the basis of a single molecular data 
alone (Ferguson, 2002). Therefore, in our study keeping the /6S , another 
faster evolving Cytb gene was used to infer the phylogenetic relationships 
among closely related and poorly known species. In amphibians, 
morphological, ecological, and/or bioacoustics characteristics have been used 
to define species (Hasan et al., 2014b). Acoustic data are important to 
understand breeding behavior of a frog as well as facilitate species tracking 
based upon its calling during field survey. The morphometric comparison of 
M. nilpahamariensis (newly described species) with M. fissipes is 
scientifically meaningful more than the comparison with M. rubra because the 
latter can be easily distinguished by its shovel like inner metatarsal tubercle 
and most importantly, it is only reported in semi-evergreen forest in 
Chittagong, SE Bangladesh (Asmat et al., 2003) rather than the northern 
region. Further, the coloration in life of any newly described species is also 
useful for identification among cryptic species. Therefore, the description of 
the new species by Howlader et al. (2015) fell short in providing these 
evidences. 


In amphibians, /6S is considered as a marker for outlining the 
taxonomic status of frog species (Vences et al., 2005), while the protein- 


coding Cytb, which is known to exhibit more rapid nucleotide substitutions 
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than rRNA genes, is considered to be more phylogenetically informative 
between conspecific populations (Koike and Matsui, 2003). Together with the 
molecular and morphology data, advertisement calls also play an important 
role in classifying cryptic species. This technique has recently been used to 
identify morphologically similar frog species (Hasan et al., 2012b; Matsui et 
al., 2014). 


Important biological and ecological information that are needed to 
accurately describe M. nilphamariensis remains lacking (e.g., natural history 
and acoustics data are necessary to know the breeding biology which support 
to define future conservation strategies) in the publication of Howlader et al. 
(2015). Notwithstanding, biological and/or ecological data are more useful 
rather than the genetic data for local taxonomists in a developing country like 
Bangladesh where genetics/molecular research is a matter of cost. Therefore, 
in this study, we analyzed the mitochondrial DNA (mtDNA) sequences of both 
16S rRNA and Cytb, and collected detailed acoustic and morphological 
information of this species. 


MATERIALS AND METHODS 


Specimen Collection. Specimens of M. nilphamariensis were 
collected from two localities in Bangladesh: Parbatipur in Dinajpur District 
and Berakhuti in Barua, Nilphamari District. 


DNA Extraction, Amplification, and Sequencing. Specimens and/or 
tissue samples were stored at the Institute for Amphibian Biology, Hiroshima 
University, the Department of Fisheries Biology and Genetics, Bangladesh 
Agricultural University, or the Museum of Herpetology Laboratory 
Bangladesh, Ichamati College, Dinajpur, Bangladesh. The methods for the 
extraction, PCR amplification, and sequencing of mtDNA 76S and Cytb 
fragments follow previous protocols (see Hasan et al., 2012a, b; Nasrin et al., 
2014). The obtained /6S and Cytb sequences are deposited in the 
DDBJ/EMBL/GenBank database (Accession Nos. LC090055—LC090057 for 
16S and LC090058—LC090063 for Cytb). The 165 and Cytb sequences of 
other closely related species were retrieved from GenBank (Table 1). 


Sequence alignment. The 16S and Cytb sequences of M. 
nilphamariensis and closely related species (Table 1) were first aligned 
separately using the ClustalW program (Thompson et al. 1994). The 
divergence of the aligned sequences (uncorrected P values) was calculated 
using MEGA Ver. 6.0 (Tamura et al., 2013) with the pairwise-deletion option 
in which all alignable sites were used for calibration, but indel sites were not 
counted. Gaps and ambiguous sites were excluded using Gblocks Ver. 0.91b 
(Castresana, 2000) with the default parameters. Gaps in the alignments were 
treated as missing data. The initial two alignments (J6S and Cytb) were 
combined into one concatenated data set, which contained a total of 1251 sites 
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(755 for 16S and 496 for Cytb), 293 of which were parsimoniously 
informative. 


Phylogenetic analyses. The phylogenetic analyses were performed 
using the maximum likelihood (ML), maximum parsimony (MP), and 
Bayesian inference (BI) methods. Nucleotide substitution models for ML and 
BI analyses were selected based on the Akaike information criterion and 
Bayesian information criterion, respectively, that are implemented in the 
Kakusan 3.0 program (Tanabe, 2007). ML analysis was performed using 
Treefinder (Jobb, 2008) and the resultant tree was evaluated by bootstrap (BS) 
analysis with 1000 replicates. BI analysis was performed using MrBayes Ver. 
3.1.2 (Ronquist and Huelsenbeck, 2003) with the following settings: number 
of Markov chain Monte Carlo (MCMC) generations = 10 million; sampling 
frequency = 100; and the first 1 million generations were discarded as burn-in. 
The number of MCMC generations and burn-in size were determined by 
checking the convergence of —log likelihood (—InL) values and tree length 
against generation number using Tracer Ver. 1.4 (Rambaut and Drummond, 
2007). Statistical support of the BI tree was evaluated by Bayesian posterior 
probability (BPP). MP was performed with 1000 bootstrap (BS) replicates 
using PAUP’ 4.0b10 (Swofford, 2003). 


Recording advertisement calls. The advertisement calls of M. 
nilphamariensis were recorded at Berakhuti, Barua in Nilphamari District (25° 
49.320' N, 88° 49.732' E) using a Sony Cybershot digital camera in video 
mode (model no: DSC-W730). Sound spectrograms were depicted using 
Avisoft-SASLab Light Software (Avisoft Bioacoustics, Germany) with 
Hamming window (111 Hz bandwidth). The calls of Microhyla species were 
long series of notes consisting of about 10 pulses (Kuramoto and Joshy, 2006). 
We measured note length, number of pulses, pulse repetition rate and 
dominant frequency. For comparison, calls of M. ornata recorded at Mudigere, 
Chikmagalur District, Karnataka State, India and those of M. fissipes recorded 
at Manchou, Pingtung County, Taiwan were analyzed. 


Morphological measurements. For morphological comparisons, the 
following 29 measurements were taken to the nearest 0.1 mm using digital 
calipers: SVL, snout-vent length; HL, head length; HW, head width; S-N, 
snout to nostril distance; N-N, inter-nostril distance; N-E, nostril to eye 
distance; ED, longitudinal eye diameter; E-E, inter-orbital distance between 
inner borders of upper eyelids; ELW, eyelid width; FLL, forelimb length; 
FHL, forearm and hand length; FAW, forearm width; HAL, hand length; F1- 
F4, lengths of Ist to 4th finger; HLL, hind limb length; FEL, femur length; 
TIL, tibia length; TFL, tarsus and foot length; FOL, foot length; T1-TS5, 
lengths of Ist to 5th toe; IMT, inner metatarsal tubercle length; and OMT, 
outer metatarsal tubercle length. Statistical analysis was performed in SPSS 
(15.0J) software (SPSS Japan Inc., Tokyo, Japan). 
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Table 1. Samples used in this study, and GenBank accession numbers. IABHU: Institute for Amphibian Biology, Hiroshiam University; 
MHLB: Museum of Herpetology Laboratory' Bangladesh; BNHS: Bombay Natural History Society; RBRL: Rondano 


Biodiversity Research Laboratory, St. Aloysius College; ZSIK: Zoological Survey of India, Kolkata. 





Accession numbers 








Species Locality Voucher no. 16S Cytb 

M. nilphamariensis Berakhuti, Barua, Nilphamari, Bangladesh IABHU 4212 LC090055 = 

M. nilphamariensis Berakhuti, Barua, Nilphamari, Bangladesh IABHU 4213 LC090056 =- 

M. nilphamariensis Berakhuti, Barua, Nilphamari, Bangladesh MHLB00206 LC090057 == 

M. nilphamariensis Parbatipur, Dinajpur, Bangladesh IABHU 22135 AB530537 = 

M. nilphamariensis Parbatipur, Dinajpur, Bangladesh IABHU 22136 AB530538 AB819032 
M. nilphamariensis Parbatipur, Dinajpur, Bangladesh IABHU 22137 AB530539 LC090058 
M. ornata Bajipe, India Released AB530627 LC090059 
M. ornata Karnoor, India BNHS 5028 AB530628 LC090060 
M. ornata Talagini, Shimoga, India RBRL 040723-04 AB530630 LC090061 
M. ornata Dharwad, Karnataka, India ZSIK-A9119 AB201188 AB201223 
M. rubra Dharwad, Karnataka, India Released AB201192 AB201224 
M. berdmorei Gombak, Malaysia IABHU21019 AB530638 AB819031 
M. cf. berdmorei Sylhet & Bandarban, Bangladesh IABHU 3864 AB530541 AB819030 
M. mymensinghensis Mymensingh, Bangladesh IABHU3899 AB543607 AB819019 
M. mukhlesuri Chittagong, Bangladesh IABHU 3880 AB543609 LC090062 
M. fissipes Huangshan, Anhui, China KUHE32943 AB201185 AB201213 
M. fissipes Thong Pha Phum, Kanchanaburi, Thailand KUHE35165 AB201186 AB201215 
M. okinavensis Ishigaki Island, Okinawa, Japan IABHU5263 AB303950 AB303950 
M. heymonsii China Living Animals AY458596 AY458596 
M. heymonsii University Malaya Campus, Malaysia IABHU 21026 AB530637 LC090063 
K. pulchra Bandarban, Bangladesh IABHU3781 AB530543 AB819033 
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RESULTS AND DISCUSSION 


Molecular phylogeny. In the phylogenetic trees that we constructed 
(Fig. 1), M. nilphamarienisis from the Dinajpur District formed a clade with 
high bootstrap support (BS: 86 for ML; 100 for MP, and 1.0 for BD, and the 
topotypic M. ornata from Bajipe, Karnoor, and Talagini in India formed 
another clade also with high bootstrap support (BS: 92 for ML; 100 for MP, 
and 1.0 for BD. The M. nilphamariensis from Dinajpur District formed a sister 
clade with the topotypic M. ornata clade with high bootstrap support (BS: 90 
for ML; 97 for MP, and 1.0 for BD. The sequence divergence between these 
two clades was 3.2% for 16S and 10.6% for Cytb. We found that all the 
specimens from India and Bangladesh formed a clade, while the Microhyla 
frogs from Southeast Asia (China, Japan, Malaysia, and Thailand) and two 
newly described Microhyla species from Bangladesh formed another clade. 
This result may reflect a sister relationship between Indian subcontinental 
species and other Southeast Asian frogs. 


The herpetofauna of Bangladesh is poorly known and identification 
of the morphologically conserved frogs of the genus Microhyla is a difficult 
task (Mahony et al., 2009; Hasan et al., 2012a). The accurate delineation of 
Microhyla species has been challenging because of their minute body size and 
high level of homoplasy as a result of the loss of pectoral girdle elements 
(Zweifel, 1986). Recently, Howlader et al. (2015) described a new species, M. 
nilphamariensis, from the northern part of Bangladesh, but Cytb data and 
acoustic information was not provided. Microhyla nilphamariensis was found 
to be close to M. ornata in both the phylogenetic and morphometric analyses. 
Because of their morphological homoplasy, a clear-cut separation among 
microhylid frogs is very difficult (Kuramoto and Joshy, 2006). In the present 
study, M. nilphamariensis formed a sister clade with the topotypic M. ornata 
from Bajipe, Karnoor, and Talagini in India, and exhibited a high degree of 
genetic divergence (3.2% for 168 and 10.6% for Cytb). It has been suggested 
that there may be more than two lineages in the currently known M. ornata 
species from India (Hasan et al., 2014b). Further, the M. nilphamariensis 
embedded within the Indian microhylid clade, which was remotely related 
with Southeast Asian microhylid frogs. Similar evolutionary relationships 
have been observed among frogs in the genus Fejervarya (Hasan et al., 
2014b). 


Generally, it is hypothesized that during the uplift of the Himalayas, 
the Bengal basin (including present-day Bangladesh) was formed between 20 
and 14 million years ago (Mya) (Alam et al., 2003; Uddin and Lundberg, 
2004). Recently, Hasan et al. (2014a) proposed that M. nilphamariensis and M. 
ornata diverged from each other approximately 19.40 (45.23-6.90) Mya. This 
event seems to have occurred before the land formation of Bangladesh and 
was continued until the late Miocene period. The common ancestor of M. 
nilphamariensis may have evolved somewhere in South Asia, and this species 
split from its near congener before 19.40 Mya (Hasan et al., 2014a). 
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Kaloula pulchra (Bandarban, Bangladesh) 
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Nucleotide substitution/site 


Figure 1: Maximum Likelihood (ML) tree based on the 1251 bp of 
mitochondrial (165 + Cytb) genes. The numbers near branches 
represent bootstrap supports for ML and MP inferences, and 
Bayesian posterior probability (ML-BPs/MP-BPs/BPP). The scale 
bar represents 0.1 nucleotide substitution/site. 


Advertisement call. Sound spectrograms and acoustic parameters of 
advertisement calls are given in Fig. 2 and Table 2. Notes of M. 
nilphamariensis were slightly longer and involved more pulses than those of 
M. ornata and M. fissipes. Pulse repetition rates were nearly the same in M. 
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Table 2: Acoustic parameters of advertisement calls of three Microhyla 
species (mean + S. D.). 


Note length Mil Puente Lower dom. Higher dom 






Species n . pulses 
(s) (pulses!) — freq, (KHz) freq, (kHz) 
M nilphamariensis y 041=001 1511+033 3613+039 16200 3.622 0.06 
M. omata M 0342005 1179+167 3167+09% 14007 1022007 
M fissipes H 036=003 1309+207 3527+150 13620% 2372001 
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Figure 2: Sound spectrograms showing advertisement call structures of M. 
nilphamariensis (A), M. ornata (B) and M. fissipes (C). 
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nilphamariensis and M. fissipes, but the range in M. nilphamariensis (35.4- 
36.5 pulses/s) did not overlap with that of M. ornata (31.6-34.4 pulses/s). 
Pulse rate of M. ornata from Assam and Meghalaya, India (Roy, 1996; 4.6 
pulses/s) and that of M. ornata from Thailand (Heyer, 1971; 53-63 pulses/s) 
differed completely from our values, and the taxonomic status of their 
materials should be examined in future studies. 


The call of M. nilphamariensis had two distinct dominant frequency 
bands at 1.65 and 3.62 kHz (Fig. 2A). Similarly, the calls of M. ornata and M. 
Jissipes exhibited two dominant frequency bands (Figs. 2B, 2C), but they were 
much lower than those of M. nilphamariensis. The ranges in M. 
nilphamariensis (1.59-1.71 and 3.52-3.70 kHz) did not overlap those of M. 
ornata (1.14-1.36 and 2.54-2.81 kHz) and M. fissipes (1.31-1.42 and 2.85-2.88 
kHz). The fundamental acoustic parameters reported for M. rubra (note 
duration 0.17 s and pulse rate 108 pulse/s) (Kanamadi et al., 1994) were 
clearly very different from those of the other Microhyla species mentioned 
above. 


Morphometric comparisons. The measurements of 29 characters of 
M. nilphamariensis, M. fissipes, and M. ornata are listed in Table 3. Microhyla 
nilphamariensis was smaller in size (SVL) = 14.74 + 2.57 mm) than topotypic 
M. ornata (SVL = 24.10 + 1.17 mm) and M. fissipes (SVL = 21.82 + 0.85 
mm). Further, M. nilphamariensis differed significantly from its close relatives 
in relation to the lengths of the first, second, and fourth fingers (Fl, F2, and 
F4), the first toe (T1), the inner metatarsal tubercle (IMT), and the outer 
metatarsal tubercle (OMT) (P < 0.04) (Table 3). 


The three species were clearly separated by principal component 
analysis (PCA) (Figs. 3A, 3B). Eigen values were 18.782 for function 1, and 
2.982 for function 2. Coefficients were large in SVL, F2, HLL, TIL, FOL, and 
TS for function 1 and in FLL and TIL for function 2. In the PCA (Fig. 3B), M. 
nilphamariensis formed a cluster that was separate from the other two species, 
while the scores of M. ornata and M. fissipes overlapped considerably (Fig. 
3A). 


Body ratios relative to SVL (e.g., head length (HL)/SVL and head 
width (HW)/SVL) and 10 other ratios of M. nilphamariensis vs. M. ornata and 
M. fissipes as well as the results of the Mann—Whitney U test between ratios 
are shown in Table 4. The results showed that M. nilphamariensis differed 
significantly from both M. ornata and M. fissipes, having smaller F1, F2, and 
T1 and larger IMT and OMT ratios relative to SVL (P < 0.04) (Table 4). 
Further, M. nilphamariensis differed from M. ornata in having smaller 
forearm width (FAW) and shorter 2nd and 4th toes (T2 and T4) relative to 
SVL, and from M. fissipes, in having smaller nostril to eye distance (N-E), 
FAW, F4, hind limb length (HLL), and tibia length (TFL) relative to SVL (P < 
0.04) (Table 4). 
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Morphologically, M. nilphamariensis was easily distinguished from 
its near congeners, having SVL of 11.7-18.3 mm, finger lengths with 1 < 2 < 4 
< 3, and rudimentary toe webbings. 





Figure 3: Results of the multivariate analyses of morphometric variability in 
three Microhyla species. Scatterplot of (A) principal component 1 
(PCA1) versus principal component 2 (PCA2) and (B) principal 
component 1 (PCA1) versus principal component 3 (PCA3) from 
the principal component analysis of M. nilphamariensis (rhombus 
[blue]), M. ornata (rectangle [purple]) and M. fissipes (triangle 
[yellow]). 


Color in life. In life, the dorsum varied from brown or slightly 
reddish with a black band-type mark running from the head to the vent. A few 
black spots were observed on the dorsal surface of the thigh and femur 
regions. The ventral surface was slightly whitish with a black mark on the 
throat and a few big black dots between the bases of the forelimbs in males of 
M. nilphamariensis (Figs. 4A, 4B). Photographs in life of M. ornata and M. 
fissipes were shown in Figs. 4C and 4D, respectively. The X-shaped mark on 
the dorsum of these three species facilitated to separate them from each other. 


Distribution. Microhyla ornata is widely distributed in poorly 
documented areas including Kashmir (Pakistan and India); Nepal; throughout 
India, including the Andaman and Nicobar Islands, Sri Lanka and Bangladesh. 
Microhyla fissipes is distributed in southern and central China (including 
Taiwan and Hainan) from southern Yunnan northward and east to Shanxi and 
Shaanxi, Thailand and Indochina and through the Malay Peninsula to 
Singapore and possibly to Bangladesh (Frost, 2015) (Fig. 5A). 


In this study, we collected individuals of M. nilphamariensis from 
Berakhuti, Barua in Nilphamari District. Berkahuti is around 7 km from Koya 
Golahut, Saidpur (also in Nilphamari District), which was considered as the 
type locality of M. nilphamariensis by Howlader et al. (2015). M. 
nilphamariensis has been found in the northern part of Bangladesh (Dinajpur 
and Nilphamari districts); however, its range may extend around the currently 
known distribution area (Fig. 5B). 
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Figure 4: 


Comparisons of three Microhyla species in life. (A) Microhyla 
nilphamariensis from Parbatipur, Dinajpur [[ABHU 22135], (B) 
Berakhuti, Barua, Nilphamari [[ABHU 4212], Bangladesh; (C) M. 
ornata from Bajipe, Mangalore, Dakshin Kannad District, 
Karnataka State, India and (D) M. fissipes from Manchou, Pingtung 
County, Taiwan. Scale bar = 10 mm. 






































Figure 5: 





(A) Map showing the approximate distribution areas of M. 
nilphamariensis, M. ornata and M. fissipes. (B) Bangladesh map 
showing the two collecting localities (open circle) and reported 
type locality by Howlader et al., 2015 (solid triangle) of M. 
nilphamariensis. 
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Natural history. We observed M. nilphamariensis either on grass in 
an open field or near shallow water of croplands. This species is abundant in 
agricultural croplands and home state forest areas. Although no other 
Microhyla species were found, many Fejervarya sp. large type were observed 
in the same locality (Berakhuti, Barua, Nilphamari in Bangladesh). 


CONCLUSIONS 


In conclusion, M. nilphamariensis can be separated from other 
microhylid frogs based on the sequences of its mtDNA (16S and Cytb), and 
acoustic and morphological data. In the future, it will be important to do more 
sampling in the whole of Bangladesh to clearly understand the distribution 
status of this frog species as well as to study the fine-scale population structure 
and genetic variation within the population. Habitat loss and fragmentation 
have reduced population sizes and have led to extinction of many amphibian 
species in human altered landscapes (Becker et al., 2007); similar 
anthropological and environmental changes are frequently occurring in 
Bangladesh (Islam et al., 2012). Therefore, the existing natural population 
might be under threat and necessary measures will need to be taken for proper 
conservation of this species. Likewise other microhylid frogs, for example, 
Kaloula pulchra and Uperodon globulosus have been marked as Vulnerable 
and Endangered, respectively, in Bangladesh (IUCN-Bangladesh, 2000). 
Therefore, the Government of Bangladesh formulated the Wildlife 
(Conservation and Security) Act, 2012 to protect threatened amphibian and 
reptiles and has declared 17 National Parks, 20 Wildlife Sanctuaries, 12 
Ecologically Critical Areas, 6 Eco-Parks, 2 Safari Parks and 2 Botanical 
Gardens that are the home of many amphibians and reptiles (Hasan et al., 
2014c). 
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Appendix: Examined specimens list 


Microhyla ornata: (13 specimens): Rondano Biodiversity Research 
Laboratory, St. Aloysius College: RBRL 040723-04, 00062617—18, 
04071113—122. 

Collection locality: Talagini and Karnoor, Karnataka, India. 


Microhyla fissipes: (15 specimens): Osaka Museum of Natural History: 
OMNH am 20028—20042. 
Collection locality: Pingtung, Taiwan. 


Microhyla nilphamariensis (9 specimens): Institute for Amphibian Biology, 
Hiroshima University: IABHU 22135—22137, 4212-4213 and MHLB 201- 
202, 205-206. 

Collection locality: Parbatipur, Dinajpur and Berakhuti, Barua, Nilphamari 
Sadar Upazila, Nilphamari, Bangladesh. 
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